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Abstract 



A modified Gross-Pitaevskii approximation was introduced recently for 
bosons in dimension d < 2 by Kolomeisky et al. (Phys. Rev. Lett. 85 1146 
(2000)). We use the density functional approach with sixth-degree interaction 
energy term in the Bose field to reproduce the stationary-frame results of 
Kolomeisky et al. for a one-dimensional Bose-Einstein system with a repulsive 
interaction. We also find a soliton solution for an attractive interaction, which 
may be boosted to a finite velocity by a Galilean transformation. The stability 
of such a soliton is discussed analytically. We provide a general treatment of 
stationary solutions in one dimension which includes the above solutions as 
special cases. This treatment leads to a variety of stationary wave solutions 
for both attractive and repulsive interactions. 



PACS: 03,75.Fi, 05.30.Jp 
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I. INTRODUCTION 



There has been a number of theoretical studies on Bose-Einstein condensation (BEC) in 
a trap that is effectively lower dimensional A three-dimensional (axially symmetric) 

anisotropic harmonic trap may be regarded as (quasi) two- or one-dimensional when uj_ <C 
u x , or u± ^> uo x respectively, at temperatures much smaller than the larger of the oscillator 
energy gaps. Recently, Burger et al. || have experimentally created dark solitons in such a 
quasi one- dimensional system of dilute Rb 87 atoms. In this experiment, u± was taken to be 
about 30 times larger than uo x . After the creation of the soliton by phase- imprinting on the 
cigar-shaped condensate, the trap is switched off, and the soliton is observed to propagate 
along the x-direction. The s-wave inter-atomic interaction between Kb 87 atoms is repulsive, 
and stable dark solitons are expected to exist in one dimension, as seen in this experiment. 
We note, however, that for Rb 85 and Li 7 , the interatomic interaction is attractive, and bright 
solitons may form. Such solitons have been studied for the quasi one-dimensional case by 
several groups @|J, using the Gross-Pitaevskii (GP) mean- field theory ||10|| . The standard 
GP theory replaces the bosonic field operator by a classical order parameter field $(r, t) 
in three dimensions, and the potential energy is given by a fourth-degree |$| 4 term in the 
energy-density functional. Perez-Garcia et al. j|] have shown how this may be reduced, for 
extreme anisotropy ^> u x ), to the quasi one-dimensional geometry with a fourth-degree 
potential energy term. It has been also shown that for the one- dimensional case with this 
attractive fourth-degree interaction potential, the solitons are dynamically stable ( for the 
repulsive case, this has been known for a long time [[Hp). One has also considered the 
limiting case for uj x — 0, when there is no confinement in the x-direction, and studied the 
free propagation of the soliton. 



A different approach has been taken recently by Kolomeisky et al. [12]. Instead of 
starting with the usual three-dimensional GP energy functional and reducing it to the quasi 
one-dimensional described above, they have derived the GP functional for truly 



lower dimensions, d < 2, using the renormalization group approach |[13|| . In particular, 
they proposed, for a repulsive interatomic interaction, a sixth-degree potential energy term 
replacing the usual fourth-degree term in the order parameter. With this modification they 
obtained a (dark) soliton solution in the repulsive case. 

In this paper, we first discuss the solutions of the modified GP formalism in one dimension 



in the stationary frame. We show that the soliton solution of Kolomeisky et al. |L2[ may 
be reproduced in the stationary frame with a repulsive interaction. Since the number of 
particles (see the definition in the next section) is generally not finite unless regulated, 
the phase cannot be obtained by a Galilean boost even though the system has Galilean 
invariance to begin with. 

Next we focus on an attractive sixth-degree interaction term in the energy functional 
for bosons. Unlike the repulsive sixth-degree interaction, the attractive interaction may not 
be obtainable as a modified GP functional in one dimension. However, the motivations for 
taking such a phenomenological attractive potential energy functional are two-fold. In three 
dimensions, using the standard fourth-degree attractive interaction energy term in the GP 
functional, it was predicted that the soliton-like solution would be stable only up to a limited 
number of bosons in the condensate B,|T4]-|T6|| , and this was verified experimentally |T7| . Of 
course, in one dimension with an attractive fourth-degree term, the soliton is variationally 
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stable as we will show in Sec. Ill D below. We want to examine the analytical solutions for 
the attractive sixth-degree interaction in the one-dimensional case. Analytical solutions for 
a more general case which includes the sixth-degree interaction have been derived earlier, 
in a different context, by Pelinovsky et al. ||18|| . They have also discussed numerically the 
stability of these solutions. Here we use the density functional approach to derive these 
analytical solutions in the stationary frame. We show that there exists a limitation on the 
particle number for stability. Since the particle number is finite in the attractive case, we 
may use the Galilean invariance to obtain the general solution in an arbitrary frame. 

It should also be mentioned that the fourth- and sixth-degree interaction (in the field 
variable) in the energy functional has been used extensively in various applications of physics 
|i"9| , with a recent one |H| examining the stability of the Li 7 condensate with an attractive 
two-body and a repulsive three-body contact interaction in three dimensions. More recently, 
the stability of an attractive BEC in the presence of a sixth-degree interaction with confine- 
ment has been analyzed by Zhang The purely sixth-degree energy functional in one 
dimension, without the confining potential considered by Kolomeisky et al. fll~2 |, is special 



because it is scale-invariant, and it has a peculiar condition on the particle number for the 
soliton in the attractive case. We further show that there are soliton as well as stationary 
wave solutions if one uses the phenomenological model for an attractive interaction. 

In Sec. II, we discuss the equations of motion in the stationary frame and their sym- 
metries. In Sec. Ill, we obtain the solutions for repulsive and attractive interactions, and 
discuss the stability of the soliton solution in the attractive case. Sec. IV provides a general 
discussion of the stationary solutions which includes all the above solutions as special cases. 



II. EQUATIONS OF MOTION 



Our starting point is Eqs. (3) and (4) of Kolomeisky et al. |12|] without a trap potential. 



The modified GP functional for a one- dimensional Bose-Einstein system is given by 

E = 7T- dx[ —— + 9 — ($$*) 3 ] , (1) 



2m J dx dx 

The first term is the kinetic energy density and the second term arises from interactions; 
here the dimensionless coupling g may be positive or negative. Note that the interaction 
term in is cubic in the density p = and not quadratic as in the standard theory. A 
rigorous derivation of this term for the repulsive case may be found in Ref. jCJ- It suffices 



to state here that in one dimension, where particles are not allowed to cross, such a term 
naturally arises when the interaction is short-range. This has been demonstrated for the 
hard-core interaction by Girardeau |22| in the thermodynamic limit, and by Lieb and Liniger 



for a repulsive delta function potential. Even with a two-body interaction which varies 
inversely as the distance between particles, the leading term in the density functional has a 
cubic dependence to the density f24}| . 

A few remarks may be made here in connection with the comment made by Bhaduri 
and Sen |25[] recently. The Lieb and Liniger analysis |23| yields the result that in the limit 
of very low density of very strong repulsive interaction, the system becomes fermionic with 
g = 1 in the functional given in Eq. (|l|). In the ground state, the density is constant and 
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the kinetic term vanishes. Hence one obtains the free fermion result for the energy per unit 
length, 

E = p 3 . 2 

2m 3 H K ' 

Obviously, in the free fermion limit no soliton solutions occur. In what follows, we keep 
both the kinetic term and the interaction term between the bosons with a general coupling 
g, and therefore are not in the free fermion limit. 

The equation of motion may be obtained from the energy functional in Eq. ([I]), and is 
given by 

ft 2 r d 2 2 _ l4 . _ , d$(x,t) 

Henceforth we shall denote the partial derivatives with respect to x by a prime and with 
respect to t by a dot. It is useful to multiply Eq. (|3|) by <&*(x, t), 

h 2 

- 7^$*$" + ®*W(x)<$> = zft$*$ , (4) 

and take the complex conjugate, 

k 2 

+ $W(x)$* = -ih$$* . (5) 

2m 

where W(x) = (h 2 /2m)gir 2 \§\ 4 . 
We now set, quite generally, 



$(x,t) = VPOM) ^ , (6) 
where p and 9 are real. It is useful to define the total particle number 

N = J dx p , (7) 

and the momentum functional 
„ — ih 



J dx [ - } = h J dx p6' . (8) 

We then obtain the equations of motion as follows. The difference between Eqs. @ and 
(||) reduces to the continuity equation for the number density p, 

p + A {p6 y = o . (9) 

Thus it is legitimate to regard the velocity of the "fluid" to be 

v(x,t) = — 9'ix.t) . (10) 
m 
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The energy density equation is obtained by adding Eqs. (f|) and (|) and dividing by 2, 

r{x,t) + W(x)p = -hp9 , (11) 
where the kinetic energy density is 

h 2 

T f x t) = ($*$" + <£>*"$) 

4m 

Note that the Eqs. (|^) and ( |TT] ) are invariant under scaling and Galilean boosts [|24 . 
Under scaling by a factor a, 

p(x,t) — > ap(ax, a t) , 

6(x,t) -> 9(ax,a 2 t) . (13) 

Then the particle number N remains invariant, while the energy and momentum transform 

as 

E a 2 E and P -> aP . (14) 
Under a Galilean boost by velocity v, 

p(x,t) — > p(x — vt,t) , 

6(x,t) -> 6(x-vt,t) + — (x - -vt) . (15) 

Then the particle number iV remains invariant while the energy and momentum transform 

as 

E ^ E + Pv + -mNv 2 , 
2 

P -> P + mNv . (16) 

We should mention a subtlety here. If we are interested in solutions involving a single 
localized object (such as the solitons discussed in the next section), then it is reasonable 
to demand that its energy and momentum should only change by a finite amount under 
a boost. From Eq. (PJ), we see that this is true only if the particle number N is finite. 
Therefore we will perform boosts only when N is finite. 

III. STATIONARY SOLUTIONS: SOME ANALYTICAL RESULTS 



We will now find some stationary solutions of Eqs. (|9]) and (|Tl|). For such solutions we 
make the following ansatz, 

$(x,t) = \J p(x) exp (-ifit/%) , (17) 

namely, we set 8(x,t) = —pt/h. Note that this ansatz is consistent with the continuity 
equation since both p and 6' are zero. The equation for p(x) is then given by 

Then the following solutions exist. 
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A. Thomas-Fermi limit 



Assuming the bulk density to be a constant, we set p(x) = po- This is an allowed solution 
with 



Po 



2mp/c/rf, (19) 



provided g is positive. This is the Thomas- Fermi result ||12|| . Computing the energy, we 
reproduce the result given in Eq. (|2|) when g = 1. 

B. Repulsive interaction 

We assume g is positive and look for solutions when the density is not a constant. The 
Kolomeisky et al. solution given in their Eq. (12) may be written in the stationary frame 
with their velocity parameter (3 = [ITJJ. The solution is given by 



/ \ k cosh(2fcc) — 1 . . 

p(x) = ) , 20 

Fy ' n^g cosh(2A;x) + 2 v ; 

with 

a = . 21 

^ 2m K ' 

Note that the density tends to po = 1 at x — > ±oo; hence it is better to define the particle 
number of the soliton to be 



; , r / m 3 f°° dt V% , 

Nsouton = / dx po - p x = - — - / — = — - In 

' ^n^fg J-oo cosht + 2 n^/g 



+ 1 
y/3-1 



(22) 



In non-linear optics this number corresponds to the total intensity of light with a constant 
background subtracted out. Since N so i iton is related to g, there exists exactly one solution 
for the system for a given value of g (up to a scaling of k). 

The above solution is called a dark soliton in optics P6| because the density goes to zero 
at x = 0, and the soliton physically corresponds to a dip in the density distribution. Note 
that the above solution is valid when g is positive, i.e., for repulsive interactions. 

One can get the energy, E of the soliton by integrating the energy density [p7 given in 



Eq. dH). Using the soliton profile given in Eq. (p0|) and using Eq. (|18D it is easy to see that 
the energy of the soliton is infinite. However if the energy of the soliton is defined as the 
difference 

E so nton = E — E , (23) 

where E is the energy of the solution with p(x) = p = k/Tc^/g, then it turns out that 
the energy of the soliton is zero. It is interesting that even in the presence of a repulsive 
interaction, it does not cost any energy to create a soliton. 



Until now we have checked the results given in Ref. [12] in the stationary frame. Next 



we consider some new solutions in the presence of an attractive interaction. 
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C. Attractive interaction 



The equation of motion for <&(x,t) in the stationary frame for p may also be solved by 
the following ansatz 

= ^_JL sech (2kx) and = - ^ . (24) 
Note that the density is normalized so that 

dx p(x) = N . (25) 
The ansatz in (p4j) satisfies the equation of motion provided that 

N 2 = - — and n = - — . (26) 
4g P 2m v ; 

Since N is real we are forced to choose g to be negative. We may therefore write 

iV 2 = — f— . (27) 
Following Ref. f26| we call this the bright soliton. We note in passing that a spatial soliton 



of the secant hyperbolic type is well-known in nonlinear optics, and has been experimentally 
realized long back pgfl , although the differential equation obeyed by the solitonic electrical 
field in optics is very different from ours. 

Since N and g are related to each other by Eq. (^7]), there exists exactly one solution 
(up to scaling) for the system given a value for either g or N. The soliton solution in this 
case is a lump, and it has its maximum at x = 0. Note the peculiar fact that the larger 
the magnitude of g, the smaller is the particle number N of the soliton. The reason for this 
is explained in the next subsection where we will show that a configuration whose particle 
number is too large is unstable to collapse. 

Interestingly, we find that the energy of the soliton in this case is zero. Since N is finite, 
we can use Eq. (|15D, to boost the solution given in Eq. (|4j). Thus 

p = sech (2k(x — vt)) , 

7T 

e = - !L + _ (I _ - vt) . ,28) 

We then find that 

E = ^mNv 2 and P = mNv . (29) 

Hence it is natural to define the mass of this soliton to be M = mN. 

We further note that there also exist stationary wave solutions apart from the above. 



Starting from the dark soliton solution given in Eq. (|20|), we can transform k ik and 



p — > —p to obtain the following stationary wave solution for negative values of g: 
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k 1 — cos(2A;a;) 

p(x) = — = ^4 , 30 

HK ' ij^g 2 + cos(2A;x) V ; 

with 

" = -^r- (31) 

Finally, note that we can apply the scaling transformation given in Eq. fll3"D for all the 
solutions given above; that would change the normalization of the density in all cases. But 
we can apply the boost transformation only for the bright soliton since that is the only case 
in which the total particle number is finite. 



D. Stability of the soliton solutions 

We will now study the stability of the solitons, particularly the bright soliton given in 
Eq. (E3), using both variational and small fluctuation analysis. Before doing that, let us 
make some preliminary comments. 

The soliton solutions obtained above for both attractive and repulsive cases are scale 
invariant. Hence the type of instability which is associated with the soliton solutions for 



certain types of nonlinear Schrodinger equations (see Barashenkov et al. [p9|l ) are not ex- 
pected here. Another form of instability which can occur is the breaking of the soliton into 
some other fundamental constituents f30fl . To check the possibility of that, let us note that 
our energy functional has a global U(l) symmetry, and hence a conserved charge associated 
with such a soliton solution. For the attractive case, the charge is given by the expression 

Q = ih fdx - = 2uN = - — . (32) 

J 2m y \g\ 

This is of course proportional to the integral of the density p(x) in our case which is the 
same in any frame. In that case the soliton solution is classically stable against the aforesaid 
complete dissociation into more fundamental constituents (mfc), because its rest energy 
which is here must be lower than mo|Q| if we assume that the mass mo of the hypothetical 
mfc must be positive. 

We will now show that in the attractive case g < 0, the bright soliton can collapse to a 
point if the particle number of a configuration is bigger than a critical number N c which is 
a function of the parameter g. We can show this as follows. Consider a trial wave function 
of the form 

= f n °, y 2 , (33) 

\7T cosh(ax) J 

where a is a free parameter (whose inverse determines the width of the configuration), and 
we set 9 = 0. The particle number for this configuration is given by N. We then find that 
the energy functional in Eq. (1) is given by 
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for this configuration. We observe that this expression does not have a minimum at any 
finite and nonzero value of a. If 



N 2 > N 2 = — , (35) 
4\g\ 



then the energy of the configuration given in Eq. ([33]) can be made arbitrarily large and 
negative by letting a — > oo; this corresponds to the density of the configuration collapsing 
to a 5-function at the origin. We thus see that there is a critical value of N c which is 
proportional to 1/ y\g\ beyond which there is no lower bound to the energy. We observe 
that the particle number of the bright soliton in Eq. (|27|) is exactly equal to this critical 
value. Since this argument is only based on a trial wave function, it does not prove anything 
about the stability of the bright soliton against a collapse. 

It is instructive to contrast the situation here with a more general form of the GP 
functional where the interaction term in Eq. ([I]) is taken to be of the form gp n instead of 
gp 3 , where we assume n > 1. If the exponent of p is taken to be n < 3 (instead of n = 3 as 
we have considered so far), and g is negative, we find from the above argument that there 
is no instability to a collapse, no matter what the values of g and N are. Namely, for any 
value of g < and N, the energy of the trial configuration given in ([53]) has a minimum at 
some finite value of the inverse width a. On the other hand, if we take the exponent n > 3 
and g < 0, the above argument indicates an instability to collapse for any value of N, no 
matter how small. Thus, an attractive interaction of the form p 3 is rather special; this is 
the only value of the exponent n for which there is a finite critical number N c , such that the 
trial configuration is stable against a collapse if N < N c and is unstable if N > N c . This 
special property of n = 3 can be traced to the fact that the coefficient of the interaction g 
is dimensionless only for this value of n. 

Let us now return to the bright soliton given in Eq. (f2~4"|), and directly study its stability 
by studying small fluctuations around it. Our main conclusion will be that there is at least 



one unstable mode. To do the analysis, let us set k = 1/2 for simplicity in ([24]). Then the 
soliton takes the form 



N 



7T cosh a; ' 

d c = P- , (36) 



8m 

with N 2 g = —3/4. We now assume small fluctuations of the form 



N B , ,s 

P = r 1" Op(x,t) , 

7T cosh a; 

fit 

= — + 56(x,t) . (37) 
8m 



The equations of motion then imply that 
dSp KN 1 d 2 sinhx d 



86 



dt urn cosh x dx 2 cosh 2 x dx 

859 nh d 2 sinh x d sinh 2 x 9 1 

— — = — — — coshx ( — - + — : — — — + k — + — — - ) op . (38) 

dt ANm dx 1 com x ox 2 cosh x 2 cosh x 2 
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To solve these equations, we eliminate 89 and define a function if)(x) as 



5p{x,t) = (coshx)~ 1/2 e - *"* . (39) 
We then obtain the eigenvalue equation 

A ift(x) = u 2 ip(x), 



where 



k 2 ,d l 1 3 x , d' 2 1 15 



A Am 2 ^dx 2 4 + 4 cosh 2 ^cfe 2 4 + 4 cosh 2 ^ ' 

The operator A is hermitian with the weight w(x) = 1; namely, for any two normalizable 
functions ^i(x) and ip2(%), 

dx w(x) 4>i(x) A 4> 2 (x) = dx w(x) (A 4>i(x)Y , (41) 



with w(x) = 1. (It is for this reason that we introduced the factor of (coshx) 1 / 2 in the 
definition of ip in Eq. (|39l) ) . Hence the eigenvalues of A, uj 2 , must be real. However, u 2 
could be negative; in that case, there would be a solution in which uj is purely imaginary 
and — iu > 0. That would be an unstable mode with the function 5p(x,t) blowing up 
exponentially as t goes to oo. We therefore need to know if A has any negative eigenvalues. 

We already know two eigenfunctions of A with uj = corresponding to the translation 
and scaling symmetries. From infinitesimal versions of these two symmetries, we find that 
the eigenfunctions of the operator A with zero eigenvalues are given by 

sinhx 

(coshx) 3 / 2 ' 

and 

1 xsinhx 
(coshx) 1 / 2 (coshx) 3 / 2 

Note that these two functions are normalizable (i.e., they would be called bound states in 
the language of quantum mechanics), and that they have one and two nodes respectively. 
Thus the operator A has quite different properties from the usual Hamiltonian operator 
H = —d 2 /dx 2 + V(x) which appears in a one- dimensional Schrodinger equation; a second 
derivative operator like H always has non-degenerate eigenvalues for bound states, and the 
corresponding eigenvalues strictly increase with the number of nodes. 

We have not been able to find an analytical expression for the eigenfunction of A corre- 
sponding to an eigenvalue u 2 < 0. However it is possible to show its existence by a variational 
argument. We use the result from quantum mechanics that for a hermitian operator A, if 
ip v (x) is a normalizable variational function, and 

= < w > 

<ip v \ip v > 
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is the corresponding variational energy, then there must be an eigenvalue of A which is less 
than or equal to E v . We now try a function of the form 

1 xsmhx 

fa[X) ~ (cOshx)V2 + 01 (coshx)^ ' ^ 

where a is a variational parameter which we assume to be real for convenience. Then the 
variational energy given in (|44 ) takes the form 



E = — + (46) 

10m 2 l + 2a + a 2 (l + vr 2 /8) ' 1 } 

We find numerically that this has a minimum value of E = -0.1350 h 2 /m 2 at a = -0.4738. 
Since this minimum value is negative, we know that there must be an eigenfunction of the 
operator A for which uo 2 < 0. This shows that the bright soliton has at least one small 
fluctuation mode which blows up as t — > oo. Since we do not know the explicit form of this 
eigenfunction, the physical interpretation of this mode remains unclear. 

Although we have found an instability classically, it is possible that the instability may 
disappear if the theory is quantized, i.e., if we impose the equal-time commutation relation 
[p(x),9(y)] = i5(x — y) |24j| . There are examples known in the literature (for instance, 
skyrmions in three dimensions) where the quantization of some collective coordinates elimi- 
nates a classical instability |HJ . However, we will not examine the question of the quantum 
stability of the bright soliton in this paper where our entire analysis has been classical. 



IV. GENERAL TREATMENT OF STATIONARY SOLUTIONS 

We will now see how the general form of the stationary solution of Eq. (§) may be found. 
We do this by mapping the one-dimensional problem considered in previous sections to that 
of a pseudo-particle problem in two dimensions. To do this let us write 

*{x,t) = e+*" l [£ 1 (x) + (47) 

Thus the density given by p = £ 2 + £| is independent of time. Eq. ([3]) then takes the form 

d 2 £i 2m/i 22 

fi + 9* P 6 , 



2 



dx 2 % 
d 2 ^2 2mn 
lb? W 



6 + gnV 6 • (48) 



These are precisely the equations of motion of a classical particle of unit mass moving in 
two dimensions with the spatial coordinates (£1,^2) an d a time coordinate x. Note that in 
the one-dimensional problem x refers to position coordinate. The total energy and pseudo- 
angular momentum of the particle are given by 

1 (d£i\ 2 1 (d&Y m/i gir 2 3 

£ = 2 fel + 2 [to) + ~W P ~ — P ' 

£ = 6^-6^- (49) 

ox ox 
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The pseudo-angular momentum L is in fact the same, up to a factor, as the momentum 
density given in Eq. @. These are conserved since the particle is moving in a central 
potential 

n&,&) = ^-^V (50) 

Since there are two conserved quantities and the particle has two degrees of freedom, the 
motion can be found exactly if the initial conditions are given. For instance, p is given in 
terms of x by the equation 

= 4{2 P S - L> - + ^V]. (51) 

We are interested in solutions of Eqs. ( Pf ) and ( |5T]) in which the particle density p 
remains finite as x — > ±00. In that case, the solution for p will generally be a periodic 
function of x; this corresponds to a stationary wave. From the form of the potential energy 
in (fi9|), we see that this will happen in the following situations. 

(i) If g > 0, the potential has a negative sixth-degree term. Hence we must have p > and 
the particle energy £ must be equal to or less than the potential barrier (this can be derived 
from Eq. (^)) in order to have p finite for all times x. Then p will always remain bounded 




by the location of the potential barrier given by po — y 2mp / gir 2 h 2 as in Eq. ([19]) • 
(ii) If g < 0, the potential has a positive sixth-degree term, so p will be finite regardless of 
whether p is positive or negative. However the possible motions are different for p > and 
p < as we will see. 

We will now discuss how the various solutions discussed in the previous section can be 
obtained as special cases from this general discussion. First of all, let us set the pseudo- 
angular momentum L = and therefore £2 = 0. Let us now consider the cases with g 
positive and negative separately. 

(i) In the repulsive case g > 0, we have already seen that we need p > and the energy 
must lie below the potential barrier so that the particle does not escape to infinitely large 
values of p. If the energy is exactly equal to the potential barrier, then we have two possible 
solutions. Either the particle stays at the top of the potential barrier for all times x (this 
corresponds to the Thomas- Fermi result in (0)), or we have a solution in which the particle 
begins at £1 = —^fpo at time x = — 00 and ends at £1 = y/po & t time x = 00. The latter 
solution corresponds to the dark soliton of Kolomeisky et al. On the other hand, if the 
energy is less than the potential barrier, then the particle will perform a periodic motion 
which repeatedly passes through £1 = 0; this corresponds to a stationary wave with the 
minimum density being zero at an infinite number of values of x. This is a new stationary 
wave solution. 

(ii) In the attractive case g < 0, there is only one kind of motion possible if p, > 0, i.e. the 
particle will perform a periodic motion passing repeatedly through £1 = 0, thus giving a 
new stationary wave. If p < 0, we have a double well potential (see Eq. (|50|), and there 
are various kinds of motion possible depending on the energy £. If £ > 0, then the particle 
again performs a periodic motion passing through £1 = 0, thereby giving the stationary wave 
described in Sec. II C. If £ — 0, then we have a bounce solution in which the particle begins 
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at £1 = at x = — oo, goes to some maximum (or minimum) value of £1 at some finite time 
and then returns to£i = 0atx = oo. This corresponds to the bright soliton described in 
Sec. II C; note that such a solution necessarily has N 2 g = —3/4 no matter what the value 
of p is. (Thus there is no bright soliton solution possible if N 2 g 7^ —3/4). If £ is negative 
but is greater than the bottom of the potential, the particle performs a periodic motion in 
which £1 always remains positive (or negative) and never reaches zero. This is a new solution 
corresponding to a stationary wave whose minimum density is nonzero. Finally, if S is equal 
to the bottom of the potential located at po = \j2mjj,/g7r 2 h 2 , then we have a solution with 
uniform density p(x) = po- 



V. COMMENTS 

We have shown that the modified one-dimensional GP mean-field theory has many so- 
lutions apart from the one outlined by Kolomeisky et al. ||12| . The soliton solution in the 



attractive case is particularly simple since the full time dependent solutions in an arbitrary 
moving frame may be obtained simply by boosting the solutions in the stationary frame. 
The question of the quantum stability of this soliton remains an open problem. 

S. G. thanks R. Rajaraman for discussions and comments. R.K.B. acknowledges financial 
support from NSERC, Canada, for this research, and the hospitality of I.M.Sc, Chennai, 
India. 
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